Entanglement renormalization in two spatial dimensions 
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We propose and test a scheme for entanglement renormalization capable of addressing large two- 
dimensional quantum lattice systems. In a translationally invariant system, the cost of simulations 
grows only as the logarithm of the lattice size; at a quantum critical point, the simulation cost 
becomes independent of the lattice size and infinite systems can be analysed. We demonstrate the 
performance of the scheme by investigating the low energy properties of the 2D quantum Ising 
model on a square lattice of linear size L = {6, 9, 18, 54, oo} with periodic boundary conditions. We 
compute the ground state and evaluate local observables and two-point correlators. We also produce 
accurate estimates of the critical magnetic field and critical exponent /?. A calculation of the energy 
gap shows that it scales as 1/L at the critical point. 

PACS numbers: 05.30.-d, 02.70.-c, 03.67.Mn, 05.50.-|-q 
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Entanglement renormalization [l| has been recently 
proposed as a real-space renormalization group (RG) 
method 0] to study extended quantum systems on a lat- 
tice. A highlight of the approach is the removal, before 
the coarse- graining step, of short-range entanglement by 
means of unitary transformations called disentanglers. 
This prevents the accumulation of short-range entangle- 
ment over successive RG transformations. Such accumu- 
lation is the reason why the density matrix renormaliza- 
tion group (DMRG) 0]- an extremely powerful technique 
for lattices in one spatial dimension - breaks down in two 
dimensions, where it can only address small systems. 

The use of disentanglers leads to a real-space RG trans- 
formation that can in principle be iterated indefinitely, 
enabling the study of very large systems in a quasi- 
exact way. This RG transformation also leads to the so- 
called multi-scale entanglement renormalization ansatz 
(MERA) [3] to describe the ground state of the system 
- or, more generally, a low energy sector of its Hilbert 
space. In a translation invariant lattice made of N sites, 
the cost of simulations grows only as log TV In the 
presence of scale invariance, this additional symmetry is 
naturally incorporated into the MERA and a very con- 
cise description, independent of the size of the lattice, is 
obtained in the infrared limit of a topological phase Q 
or at a quantum critical point [l|, 0, 0, H, S [IS EH ■ 

While the basic principles of entanglement renormal- 
ization are the same in any number of spatial dimensions, 
most available calculations refer to ID models. Numer- 
ical work with 2D lattices incurs a much larger compu- 
tational cost and has so far been limited to exploratory 
studies of free fermions Q and free bosons [sj and of the 
Ising model in a square lattice of small linear size L < 8 
[l^ . It must be emphasized, however, that the approach 
of Refs. 0, @] relies on the gaussian character of free 
particles and can not be generalised to the interacting 
whereas the results of Ref. [l2| were obtained by 



In this paper we present an implementation of the 
MERA that allows us to consider, with modest computa- 
tional resources, 2D systems of arbitrary size, including 
infinite systems. In this way we demonstrate the scalabil- 
ity of entanglement renormalization in two spatial dimen- 
sions and decisively contribute to establishing the MERA 
as a competitive approach to systematically address 2D 
lattice models. The key of the present scheme is a care- 
fully planned organization of the tensors in the MERA, 
leading to simulation costs that grow as 0{x^^), where x 
is the dimension of the vector space of an effective site. 
This is drastically smaller than the cost 0(x^^) of the 
best previous scheme 0, @, • We also demonstrate the 
performance of the scheme by analysing the 2D quan- 
tum Ising model, for which we obtain accurate estimates 
of the ground state energy and magnetizations, as well 
as two-point correlators (shown to scale polynomially at 
criticality), the energy gap, and the critical magnetic field 
and beta exponent. Finally, we discuss how the use of dis- 
entanglers affects the simulation costs, by comparing the 
MERA with a tree tensor network (TTN) 13]. 

2D MERA. — Let us consider a square lattice Co 
made oi N = Lx L sites, each one described by a Hilbert 
space V of finite dimension d. The proposed 2D MERA 
is characterized by the coarse-graining transformation of 
Fig. ([1]), where blocks of 3 x 3 sites of lattice Cq are 
mapped onto single sites of a coarser lattice Ci. This is 
achieved in three steps: first disentanglers u are applied 
on the four sites located at the corners of four adjacent 
blocks; then disentanglers v are applied at the boundary 
between two adjacent blocks, transforming four sites into 
two; finally, isometrics w are used to map a block into a 
single effective site. In this way, tensors u,v and w [14] , 
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exploiting a significant reduction in computational cost 
that occurs only for small 2D lattices. 



transform the state |^'o) G V®" of the lattice Cq in which 
we are interested (typically the ground state of a local 
Hamiltonian Ho) into a state l^-i) e V®^/^ of the effec- 
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FIG. 1: (Color online) Entanglement renormalization scheme 
for a square lattice. A block of 3 x 3 sites of lattice £t-i (i) 
is mapped onto one site of Lt (v). The RG transformation 
involves (ii) applying disentanglers u between the corners of 
adjacent blocks followed by (iii) disentanglers v which act 
across the sides of adjacent blocks and (iv) isometrics w which 
act within a block. Tensors u, v and w have a varying number 
of incoming and outgoing indices (vi) according to Eq. [T] 



tive lattice £i through the sequence 

1*0) A |vi,^) A 1*;,'} A 1*1). (2) 

To understand the role of these tensors, it is useful to 
think of the state j^o) as possessing three different kinds 
of entanglement: short-range entanglement residing at 
the corners of four adjacent blocks, short-range entangle- 
ment residing near the boundary shared by two blocks, 
and long-range entanglement. Then the disentanglers u 
and V are used to reduce the amount of short-range en- 
tanglement residing near the corners and boundaries of 
the blocks. In other words, in states |*o) and |*q) in- 
creasing amounts of short-range entanglement from j^o) 
have been removed. This fact facilitates significantly the 
job of the isometry w, namely to compress into an ef- 
fective site of £i those degrees of freedom in a block 
that still remain entangled (now mostly through long- 
range entanglement) with degrees of freedom outside the 
block. Thus, the resulting state |*i) still contains the 
long-range entanglement of |*o)j but most of its short- 
range entanglement is gone. We complete the above con- 
struction by noticing that a d-dimensional space V is of- 
ten too small to accommodate all the relevant degrees of 
freedom left on a block. Accordingly, we shall describe 
the effective sites of L\ with a space of larger dimension 
X- This dimension x determines both the accuracy and 
cost of the simulations. 
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FIG. 2: (Color online) Spontaneous and transverse magneti- 
zations (cTa;) and (cTz) as a function of the applied magnetic 
field A and for different lattice sizes L. Results for small sys- 
tems correspond to exact diagonalization whilst results for 
larger systems were obtained with a x = 6 MERA. As L in- 
creases, the magnetizations are seen to converge toward their 
thermodynamic limit values. Results for L = 54 could not be 
visually distinguished from results for L = 18 and have been 
omitted in the plot. As it is characteristic of a second order 
phase transition, for large L both magnetizations develop a 
discontinuity in their derivative, with {p^) (the order param- 
eter) suddenly dropping to zero at the quantum critical point 
(see Fig. El. 



The transformation of Fig. [T] can now be applied to 
lattice £i, producing a coarser lattice L^- More gener- 
ally, if £o is finite, O(logiV) iterations will produce a 
sequence of lattices {£oj -Ci, £2, ■ ■ • , -Ctop} where the top 
lattice /^top contains only a small number of sites and 
can be addressed with exact numerical techniques. Thus, 
given a Hamiltonian Hq on £0, we can use the above 
RG transformation to obtain a sequence of Hamiltonians 
{-ffo, Hi,H2, ■ ■ ■ , -fftop}, then diagonalize iJtop to find its 
ground state |*top), and finally recover the ground state 
1*0) of Hq by reversing all the RG transformations: 

|*top)-----|*2)-|*l)-|*0). (3) 

This is precisely how the MERA is defined. Specifically, 
the MERA for |*o) is a tensor network containing (i) a 
top tensor, that describes |*)t„p, and (n) O(logA^) lay- 
ers of tensors (disentanglers and isometrics), where each 
layer is used to invert one step of the coarse-graining 
transformation of Fig. [1] according to the sequence ^ . 

The technical details on how to numerically optimize 
the disentanglers and isometrics of the MERA to approx- 
imate the ground state |*o) of Hq are analogous to those 
discussed in Ref. Q for a ID lattice and will not be re- 
peated here. Instead, we focus on the key aspect that 
makes the present 2D scheme much more efficient than 
that of Refs. 0, For this purpose, we consider an 

operator Oq whose support is contained within a block 
of 2 X 2 sites of lattice Cq. Direct inspection shows that, 
no matter where this block is placed with respect to the 
disentanglers and isometrics of Fig. [U the support of 
the resulting coarse-grained operator Oi is also contained 
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FIG. 3: (Color online) Magnetizations (p-^) and (cr^) as a 
function of the applied magnetic field A for different values 
of the refinement parameter x- Left: Spontaneous magne- 
tization (cTi,) for L = 54. Data fits of the form {g^) ~ 
(A — Xc)^" near the critical point give a critical magnetic 
field Ac = {3.13,3.09,3.075} and critical exponent /3c = 
{0.320,0.321,0.323} for x = {2,4,6}. Current Monte Carlo 
estimates are Ac = 3.044 and fic — 0.326 [l^. Thus accuracy 
increases with x- Right: Transverse magnetization (ctz) for 
L = a. TTN results for large x ^re taken as the exact solution 
(see Fig. [S}. Whilst a x = 2 MERA produces significantly 
different values, results for x = 3 are already very similar and 
those for X = 6 MERA agree with the TTN solution on at 
least 3 significant digits. 



within a block of 2 x 2 sites of £i, and the same holds for 
any subsequent coarse-graining. This is in sharp contrast 
with the 2D scheme of Refs. [3, where the mini- 

mal stable support of local observables (or 'width' of past 
causal cones) corresponded to blocks of 3 x 3 sites. In 
the present case, much smaller objects (operators acting 
on 4 sites instead of 9 sites) are manipulated during the 
calculations, resulting in the announced dramatic drop 
in simulation costs. 

Benchmark calculations. — We have tested the pro- 
posed scheme by investigating low energy properties of 
the quantum Ising model with transverse magnetic field, 



{r.r') r 



(4) 



on a square lattice with periodic boundary conditions (lo- 
cal dimension d = 2). First of all, we consider a sequence 
of lattices with increasing linear size L = {6,9,18,54}. 
For each of them, a MERA approximation to the ground 
state of i?i3i„g for different values A £ [0, 5] of the trans- 
verse magnetic field is obtained using % = 6. Computing 
the ground state for L — 54 and critical transverse mag- 
netic field takes ~ 4 days on a 3GIIz dual-core desktop 
PC with 8Gb RAM when starting from a randomly ini- 
tialized MERA 15|. Fig. [2] displays the expected value 
of the parallel and transverse magnetizations, both of 
which show characteristic signs of a second order phase 
transition as L increases. We emphasize that since the 
simulation costs grow only as the logarithm of L, it is 
straightforward to increase the system size until e.g. fi- 
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FIG. 4: (Color online) Top: The energy gap as a function of 
the transverse magnetic field A, computed by exact diagonal- 
ization for small system sizes L = {2, 3, 4} and with a x = 6 
MERA for L — {6, 9}. The gap scales as 1/L at the critical 

magnetic field. Bottom: Two-point correlators {a^^a^ ')c at 
criticality and for different values of x- The scale invariant 
MERA produces correlators that decay polynomially with the 
distance s = \r — r'\. As x increases their asymptotic scaling 
approaches l/s^^"^ with rj = 0.03±0.01 [l3|. Correlators have 
been computed at distances 5 = 3* for fc = 0,l,2,..., where 
they can be evaluated with cost 0(x^^). For comparison, we 
have included correlators obtained with a D = 2 and D = 3 
iPEPS i_13]. The latter are very accurate for s = 1,2 but 
decay exponentially after a few sites. 



nite size effects become negligible on local observables. 

Fig. [3] shows how the parallel and transverse magne- 
tizations change with increasing X: for L = 54. Since 
the cost of the simulations grows as 0{x^^), only small 
values of x can be considered in practice. However, with 
X = 6 one already obtains estimates for the location of 
the critical point and the critical exponent /3 that already 
fall within 1% of the best Monte Carlo results [l6| . 

By using the MERA to represent a two-dimensional 
subspace and minimizing the expectation value of i?isi„g, 
we obtain the system's energy gap AE. Fig. 2] shows 
AE as a function of the transverse magnetic field and 
system size. Notice that at the critical point the gap 
closes with the system size as 1/L (dynamic exponent z — 
1). Two-point correlators can also be extracted. Fig. [4] 

shows the correlator {a'^x^a^x '}c = (clfVi'" ') — (crlT' ) (crlT ') 
along a row or column of the lattice, obtained using the 
scale invariant algorithm (lOj . which directly addresses 
an infinite lattice at the critical point. 

Role of disentanglers. — In order to highlight the 
importance of disentanglers, we have also performed sim- 
ulations with a tree tensor network (TTN). This corre- 
sponds to a more orthodox real-space RG approach where 
the block of 3 x 3 sites in Fig. [T]is directly mapped into 
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renormalization scheme for the square lattice and demon- 
strated its scalabiUty by addressing the quantum Ising 
model on systems of linear size L — {6,9, 18,54}, with 
cost 0(x^^ logi), and on an infinite system at criticality, 
with cost 0{x^^)- The key of the present approach is 
the use of two types of disentanglers that remove short- 
range entanglement residing near the corners and near 
the boundaries of the blocks while leading to narrow 
causal cones of 2 x 2 sites. Similar schemes can be built 
e.g. for triangular, hexagonal and Kagome lattices 19]. 

The authors thank Roman Orus, Luca Tagliacozzo and 
Philippe Corboz for comments. Support from the Aus- 
trahan Research Council (APA, FF0668731, DP0878830) 
is acknowledged. 



FIG. 5: (Color online) Energy error as a function of the refine- 
ment parameter x for finite systems of different sizes and for 
infinite systems. In absence of an exact solution for ground 
state energies, the errors are defined relative to the results 
obtained with (i) av = 60 TTN, (ii) a x = 9 MERA, (ui,iv) 
a D = 3 iPEPS For finite systems (i,ii), the MERA is 

compared against the TTN. The double a;-axes for Xmbra and 
Xttn have been adjusted so that they roughly correspond to 
the same computational cost. For L = & the TTN is more ef- 
ficient whilst for L = 9 the MERA already gives significantly 
better results. Comparison between MERA and iPEPS re- 
sults for (iii) an infinite system off criticality and (iv) an infi- 
nite system at criticality shows very similar accuracy between 
X = 3 MERA and D = 2 iPEPS, whereas D = 3 iPEPS gives 
a lower (better) energy than x = 6 MERA. 



an effective site without the use of disentanglers. Recall 
that a 2D ground state typically displays a boundary law, 
Si ^ I, for the entanglement entropy Si of a block of Z x Z 
sites. To reproduce this boundary law with a TTN, one 
needs to increase the dimension x each step of the 
coarse-graining. Specifically, Xttn must grow doubly ex- 
ponentially with the linear size L of the lattice. On the 
other hand, the cost of manipulating a 2D TTN grows 
only as a small power of Xttn- As a result, much larger 
values of x can be used with a TTN, leading to a very 
competitive approach for small lattice sizes [ISj. Fig 0(1 
and ii) compares the performance of the MERA and the 
TTN in lattices of size 6x6 and 9x9. It shows that a 
TTN is more efficient than the MERA in computing the 
ground state of the 6x6 lattice; however, this trend is 
already reversed in the 9x9 lattice, where the cumulative 
benefit of using disentanglers clearly outweighs the large 
cost they incur. Disentanglers, by acting on the bound- 
ary of a block, readily reproduce the entropic boundary 
law (for any value of x) a-nd allow us to consider arbi- 
trarily large systems. Fig. \5\ (iii and iv) shows results for 
an infinite lattice near and at criticality. 

To summarize, we have proposed an entanglement 
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The tensors of the MERA are called disentanglers it or 
isometries w depending on whether they are in charge 
of eliminating short-range entanglement or of mapping a 
block of sites into a single site [l|,[l, [Hj • This distinction is 
somewhat arbitrary: one can consider tensors that fulfill 
the two roles simultaneously, such as tensor v in Fig. [T] 
that we still call disentangler. All these tensors must be 



isometric, that is u'u — I, 



= I. The her- 



mitian conjugation (|) in Eq. (ITJ appears for consistency 
with previous references. 

Calculations for x = 6 are achieved by using a disen- 
tangler u with X = 4 on selected indices. The computa- 
tion time is reduced to a few hours per point by re-using 
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